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ABSTRACT IR 


This study is an extension of previous statistically 
oriented research at the Naval Postgraduate School. The 
method of Model Output Statistics is used to predict open- 
ocean visibility employing stepwise-selection, multiple 
linear regression. The visibility predictand is specified 
Categorically with comparisons made to a previous probabil- 
istic approach.  Predictors include direct and derived 
model output parameters provided by the U.S. Navy's Fleet 
Numerical Oceanography Center (FNOC), Monterey, California. 
EEcSt 138,000 North Pacific Ocean (30°-60°N) synoptic ship 
reports at 0000 GMT from June 1976 and 1977, July 1979, 
and August 1979 were used as both dependent and independent 
data sets. Visibility equations for both analysis-time 
and 24- and 48-hr prognostic times are developed, and are 
verified uSing percent correct, Heidke skill score, and 
bias. Levels of skill are less than desirable for opera- 
tional use. Important predictor parameters are found to 
be sensible and evaporative heat fluxes, meridional wind 
component, sea-level pressure, air/sea temperature differ- 
ence, relative humidity, an FNOC fog probability parameter 
and a visibility parameter derived from a marine aerosol 
model. Other experiments concerning weighted least squares 
predictand transformations and R* deflation are briefly 


described. 
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I. INTRODUCTION AND BACKGROUND 


Visibility is an important meteorological variable that 
can have a significant impact on the safety of maritime 
operations. Naval activities such as amphibious assault, 
underway replenishment and air operations can be greatly 
restricted under conditions of low visibility. Civilian 
operations can suffer also. In most cases poor visibility 
at sea is due to the occurrence of fog. The economic, mili- 
tary and human losses associated with United States Naval 
Operations attributable to fog are well documented by Wheeler 
and Leipper (1974). Thus accurate forecasts of fog, or more 
generally, marine visibility, would be of great benefit to 
the military and civilian communities. 

Earlier research into this problem at the Naval Post- 
graduate School (NPS), Monterey, California, using statistical 
methods, was conducted by Van Orman and Renard (1977), Quinn 
(1978), and Ouzts and Renard (1979), who all applied regression 
techniques to forecast the occurrence of fog with some degree 
of skill. Research into forecasting visibility directly, but 
using a very limited set of parameters and data, was conducted 
by Schramm (1966). Further work by Nelson (1972) used a 
larger data set and investigated new parameters. More recently 
the work by Aldinger (1979) continued research into determining 
those parameters which are statistically correlated with marine 


visibility. In addition, using a probabilistic approach, 
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Aldinger derived analysis-time linear regression equations 
which show a reasonable degree of probabilistic skill. He 
also expanded the evaluation of these equations to categori- 
cal estimates using Threat Score, Heidke Skill Score and 
percent correct. In addition, he adapted a scoring awards 
matrix to the verification which enhances the skill by giving 
partial credit to forecasts that are close to the observed 
ESEGgory. 

This study continues the statistical regression work on 
visibility analysis/forecasting, but uses a categorical 
approach rather than a probabilistic one. New predictor 
parameters are investigated and prognostic, as well as 
analysis-time, equations are derived. In addition, more 
attention is given to interpreting the statistical methods 


used. 


LI 





TAR OBJECTIVES 


The primary objective of this study was to expand on 
previous NPS visibility research using numerical-model output 
parameters from the Fleet Numerical Oceanography Center 


(FNOC) ‚I 


Monterey, California to diagnose and predict marine 
visibility over the open ocean by statistical means. The 
method of model output statistics (MOS) (see Glahn and Lowry, 
1972) was used to predict visibility categories directly as 
opposed to using a probabilistic approach. 

Within the primary objective, more specific goals to be 
achieved were to: 

(1) Develop statistical diagnostic (analysis-time, or Tau 
Y hr) and prognostic (forecast-time, or Tau 24 hr, 48 hr) 
visibility equations using stepwise multiple linear regression; 

(2) test several types of categorical schemes; 

E rtestovarrious forms of the visibility predictand 
in the regression program; 

(4) test predictor parameters not previously used in NPS 
visibility research; 

(5) compare the categorical approach to the probabilistic 
approach as used by Aldinger (1979); 

(6) test methods of regression other than the least- 


Squares linear type. 


lrormerly called the "Fleet Numerical Weather Central". 
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III. DATA 


A. AREA 

The area of study was limited to a region of the North 
Pacific Ocean located approximately between 30° and 60°N and 
from 145°E to 130°w. The actual area was restricted in size 
from the limits mentioned in order to reduce the number of 
land-influenced grid points used in computing derivatives 
applicable at marine grid locations. Also, this was done to 
eliminate, as much as possible, any orographic influences on 
visibility. The study area is shown in Figure 1 on a polar 
Stereographic projection, the grid points of which correspond 
to those of the standard FNOC 63x63 grid (with a mesh size 
of 381 km at 60°N). The entire FNOC grid is shown in Figure 2 
with an outlined area from which FNOC's model output parameters 
were extracted. This study area is the same as that used 
for recent statistical studies of marine fog and visibility 


at NPS. 


See SELECTION OF TIME PERIOD 

Data from the months of June, July and August only were 
used in this study. The frequency of fog - (and thus visibility) 
related maritime casualties reaches a peak during the Northern 
Hemisphere summer months (Figure 3). Therefore, this period 


is one of primary operational significance. 
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Study area on polar stereographic projection. 


Figure 1. 








Figure 2. Fleet Numerical Oceanography Center's 63x63 


ETI vith outline of North Pacific Ocean 
rectangular grid area used in study. 
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Only 0000 GMT synoptic ship report data were used as 
this ensured that daylight was present throughout the study 
area, thus allowing more accurate visibility observations 
than if nighttime observations were included. 

Model output parameter data from FNOC were taken from 
0000 GMT for use in analysis-time equations. However, in 
prognostic equations 1200 GMT parameters also were used. 

Diagnostic (Tau 0 hr) equations were developed from 
combined June 1976 and June 1977 data using analysis-time 
data only. In addition, equations for Tau 0, 24 and 48 
hrs were developed from July 1979 data using both analysis- 


time and prognostic-time parameters. 


E SYNOPTIC WEATHER REPORTS 

The synoptic weather reports used in this study were 
provided by the Naval Oceanography Command Detachment cow 
located with the National Climatic Center at Asheville, North 
Carolina. 

The total number of observations available in the area 


of Figure l is as follows: 


June 1976 (Tau O0) 4217 
June 1977 (Tau 0) 5044 
July 1979 (Tau 0) 4079 
(Tau 24) 4095 
(Tau 48) 4102 


“Formerly called the "Naval Weather Service Detachment”. 
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August 1979 (Tau 0) 4727 
(Tau 24) 4520 
(Tau 48) 4421 
The actual number of cases varied slightly from the numbers 
given above depending on experiments being performed. 
All synoptic reports from the June data sets were put 
through a quality control check by Aldinger (1979) to 
ensure a certain degree of compatability among present weather 
and visibility codes, in conformance with the Federal Meteoro- 
eal Handbook No. 2 (U.S. Depts. of Commerce, Defense, 
and Transportation, 1969). All data sets including July and 
August 1979 data were quality-control checked by the National 


Climatic Center, Asheville, N.C. 


D. INTERPOLATION SCHEME 

All model output parameters, whose positions are within 
the FNOC grid, were interpolated to the ship positions from 
which the synoptic observations were obtained. The interpo- 
lation method used is a natural bicubic spline curvilinear 
scheme. This scheme and its documentation are available at 
the NPS W.R. Church Computer Center where all the computer 


computations for this study were accomplished. 


E. PREDICTOR PARAMETERS 
l. Model Output Parameters (MOP's) 
A total of 22 analysis- and prognostic-model parameters 


were provided by FNOC. They were generated from the Mass 
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Structure Analysis model, the Primitive Equation (P.E.) 
model, and the Marine Wind model [U.S. Naval Weather Service, 
1975]. In addition, 79 other parameters were developed from 
the original set. Brief descriptions of all of these 
parameters are listed in Appendix A. 
2. Climatology Parameter 
The only climatology factor used as a parameter in 
this study is the fog climatology developed by the National 
Climatic Center [Guttman, 1978]. A suitable visibility clima- 
tology was not available at the time of this study. 
3. Interactive and Modified Parameters 
Interactive parameters were formed in this study by 
using the product of two different parameters. They have 
been used to account for possible physical interactions between 
variables. Other parameters, Called "modified", are simply 
the square, or the square root, of an MOP. A decision as to 
which variables to combine or modify out of an almost un- 
limited number of possibilities is a difficult task. There- 
fore, four of the parameters chosen here were taken from a 
previous study by Ouzts (1979). The remainder were chosen 
by combining or modifying those parameters which contributed 
significantly to explaining the variance of the predictand, 
in one or more experiments of this study. 
4. Binary Parameters 
This type of parameter is commonly used by the 


Techniques Development Laboratory of the National Weather 
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Service, Silver Springs, Maryland. A binary parameter 
is formed from an MOP by choosing one or more critical values 
of that MOP which, when equaled or exceeded, gives the binary 
a value of one; otherwise the binary has a value of zero. 
Here again, a seemingly infinite number of parameters is 
possible, but the set of binary parameters was limited to 
14 in this study. 
5. Beta Visibility Parameter 

The information for the computation of this parameter 
was supplied by Dr. A. Goroch? of the Naval Environmental 
Prediction Research Facility. The computation uses a marine 
aerosol model developed for the United States Navy to test 
electro-optical system performance. 

Apparently no formal documentation is available on 
the development of this model. However, Nounkester (1980) 
refers to this model and states that it was developed by 
modifying an empirical model proposed by Wells, et al., (1977). 
The modifications were made by B. Katz of the Naval Surface 
Weapons Center, White Oak, Maryland; L. Ruhnke of the Naval 
EL--orch Laboratory, Washington, D.C.; and M. Munn of the 
Lockheed Research Laboratory, Palo Alto, California. 

The aerosol model computes extinction coefficients and 
ranges at various wavelengths, as affected by molecular 


Scattering and absorption, aerosol extinction and weather. 


3 
Personal communication. 
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Only the visual range was of interest in this study, so 
only that portion of the model was used. 

As input, the FNOC model output surface windspeed and 
relative humidity, and present weather code were supplied. 
Then, a parameterized visibility was computed, herein called 
beta visibility (BVIS). Since two relative humidity parameters 
were available, RHR and RHX, two beta visibility parameters 
could be computed, BVISR and BVISX. 

Because the present weather code was not available 
at prognostic times, beta visibility could not be computed 
at tau 24 and tau 48. However, since the aerosol extinction 
itself was expected to correlate well with observed visi- 
bility, a modified beta visibility parameter was formed by 
simply omitting the weather code input. This modified beta 
visibility (MBVIS) could then be used at prognostic times. 
The method produced a less accurate parameter, but one that 
still correlated well with observed visibility. The methods 
used for computing the BVIS and MBVIS parameters are given 


in Appendix B.3. 
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IV. PROCEDURE 


A. REGRESSION SCHEME 

A computer program for stepwise multiple linear regression 
using the method of least squares was used to derive the 
Visibility equations. The program used is one of the UCLA 
BMDP series, namely BMDP2R [UCLA, 1979]. 

In this program the dependent variable (predictand) is 
specified, then independent variables (predictors) are entered 
(forward stepping) or removed (backward stepping) based on a 
statistical F-test with given F-to-Enter (4.0) and F-to- 
remove (3.9). The first predictor selected in forward stepping 
is the predictor variable with the highest F-to-enter.  Suc- 
ceeding steps enter variables in the same manner. At each 
Step the variables already entered into the equation are 
reevaluated and could be removed by backward stepping if they 
fail to exceed the minimum F-to-remove value. 

If a variable being considered for entry reflects a strong 
linear combination with any of the variables already entered, 
it may cause computational difficulties, and the BMDP2R 
program will reject it if its tolerance value equals or 
exceeds 0.01. The program continues stepping until all 
variables are used, or until no further variables meet the 
F-to-enter value. A further definition of the statistics 


used is included in Appendix C. 
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Another regression routine available is BMDP9R, called 
All Possible Subsets Regression. Rather than performing 
a screening regression as in BMDP2R this program considers 
all possible combinations of predictor variables to achieve 
the highest possible Be value (explained variance). This 
program was used for a few experiments. Some of the com- 
puted subsets did manage to attain a higher Re value than 
that achieved by screening regression, but these R^ values 
were only marginally higher and have doubtful significance. 
Thus, the results achieved by this method did not justify 


the excessive computer time involved, and so it was abandoned. 


B. CATEGORICAL APPROACH 

Previously at NPS, Aldinger (1979) developed analysis- 
time visibility regression equations based on a probability 
approach. Equations were developed to estimate the probability 
of occurrence of each of several visibility code groupings. 
In this study a categorical approach was used. Several schemes 
for grouping visibility codes into different categories were 
used. In order to have a visibility value for the predictand 
the midpoint value of the visibility range for each observed 
category was used. For example, if a category included synop- 
tic codes 90-93 the visibility range would be 0-1 km, and the 
visibility predictand was assigned the value of 0.5 km. An 
exception to this rule was made for the highest visibility 


category. Since this category has no upper limit, several 
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arbitrary visibility values were assigned to the predictand 
depending on the categorical scheme involved. A list of 
the synoptic visibility codes used to determine the 
visibility categories can be found in the Federal Meteor- 
ological Handbook No. 2 [U.S. Depts. of Commerce, Defense 
and Transportation]. 

The regression equations so developed yield continuous 
visibility values (in kilometers) which can be used 
directly, or perhaps more appropriately, can be used to 
specify the selected category. The latter method is used 
in this study for verification purposes. 

Since there are only ten reported synoptic visibility 
codes, with each code representing a range of visibility, 
the maximum number of defined categories is limited to ten. 
Using the maximum number of categories allows the greatest 
visibility æsolution. However, there is some inaccuracy 
involved in visibility reporting that is related ce EE GB} 
server's ability to discriminate between different visibility 
ranges. Therefore, categorical schemes were developed which 
combined several observed codes into one category. This 
approach provides a wider visibility range for each category 
and partly compensates for observer error. It is reasoned 
that an observer should be able to distinguish between a few 
larger visibility ranges better than a larger number of smaller 
visibility ranges. Of course, with fewer categories some 


visibility resolution is lost. In the extreme case, a scheme 
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with only one category, which includes all visibility values, 
would not be affected by observer error, and all regression 
estimates would be perfect. However, such a scheme obviously 
would be useless. Therefore, some tradeoff between accuracy 
and resolution should be made. In this study schemes involving 
five and ten categories were tested. 

Tau 0 equations were developed for all categorical schemes 
meemecomoined June 1976 and June 1977 data. The predictor 
parameters considered in the equations are listed in Appendix 
Apart 1. 

Analysis-time (Tau = 0 hr) and prognostic (Tau = 24 and 
48 hr) equations were developed from July 1979 data. Prog- 
nostic equations at 24 hr and 48 hr only were developed so 
that the verification times would correspond to 0000 GMT. 
However, MOP's from 00, 12, 24, 36, and 48 hr were used. The 
parameter list used to develop these equations is located in 


Appendix A, part 2. 


C. EQUATION TRUNCATION AND VERIFICATION 

The BMDP2R regression routine enters a new variable at 
each step, increasing the ES value each time, thus fitting 
the equation better to the dependent data. After a certain 
number of steps, however, the incremental increase in R? per 
step may have little or no significance when the equation is 
applied to independent data. For this reason it was decided 


to truncate each equation before entering a variable which 
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does not increase the E value by a rounded value of 1$. 
In general this produced an equation with four to six varia- 
bles. More will be said on this topic later. 

Two scoring methods were used to describe the skill of 
each final regression equation. These two methods consist 
of computing the percentage of correct forecasts and Heidke 
Skill score for each equation. The formula for computing these 
scores is given in Appendix D. The continuous visibility 
output from a regression equation lies within the visibility 
range of a particular category. This particular category is 
considered to be the one estimated by the regression equation. 
The number of times each category is thus estimated is com- 
pared to the number of observations of each category for 
Scoring purposes. 

All equations were verified against the dependent data 
from which they were derived. In addition, all five-category 
equations were verified against independent data. Equations 
developed from combined June 1976 and June 1977 were indepen- 
dently verified using July 1979 data, and equations developed 
from July 1979 data were verified using August 1979 data. 
Unfortunately, the lack of availability of MOP fields and 
observational data prevented the independent verification of 
June equations with other June data, and July equations with 
other July data. 

Another scoring technique applies a scoring matrix 


developed by Aldinger (1979) and applied to the five-category 
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scheme. The matrix applies weights to the number of esti- 
mates of each category in order to give some credit for 
nearly correct estimates. This matrix, called the NPS awards 
matrix, is further described in Section V.C.3. 

In addition, a distribution measure, called bias, is 
nated for each category. Bias represents the ratio of 
the number of forecasts to the number of observations of each 


category. 
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A. CATE 


T: 


as defin 
Category 
Number 

I 

II 

TTI 


IV 


VI 
VII 
EXIIT 


IX 


V. EXPERIMENTS, RESULTS, 


GORICAL SCHEMES 


Tenc C acegJgorFy Scheme: LOCATA 


This scheme uses ten categories of the predictand 


ed below. 


Observed 


Visibility Code 


90 
351 
52 
93 
94 
95 
96 
97 
98 


99 


Visibility 
Range (km) 

= 005 
OOS co <€ 0. 
0.2 EO <O. 
OST EOS. 
0 tons 2 
ZU EOC. 
4.0 to <10. 
POTO EO 20. 
20:0 EO 50. 

>50.0 


DISCUSSION 


0 


UN 


LO 


35: 


Ta 


Value of 
Predictand (km) 


2025 


125 


39 


29/52 


A Tau 0 equation was developed from combined June 


1976 and June 1977 data and verified on the dependent data. 


All values, except for regression coefficients are given to 


two deci 


mal places. 
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Coefficient Predictor 
-354.558 
t 1.346 EHE 
+ 0.388 BVISR 
+ 0.358 PS 
+ 5.174 SEHE 
+ 1.380 ASTDR 
- 2.938 VCMP1 
R^ - .25 
Dependent Verification: Percent Correct = 40 
Skill Score zb 
Category I Ed ITI IV V VI VII VITI IX X 
Bias USO Ol ‚ori re 56 1.60 1.46 ٠.01 


The scores for this scheme are relatively low. The 
bias values indicate that the highest category and the lowest 
six categories are observed far more often than selected by 
the regression equation. On the other hand, categories VIII 
and IX were selected much more often than they were observed. 

2.  Ten-Category Scheme:  l10CATB 

It was felt that the arbitrarily selected midpoint 
value of 75.0 km for category X in 10CATA was too high, 
thus causing a poor fit of data in the regression equation. 


Therefore, this category was changed in 10CATB, as follows. 
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Category Observed Vesıbılity Value of 


Number vireibility Code Range (km) Predictand (km) 
X 99 550 50 


All other categories, I through IX, were defined the 
same as in LOCATA. The Tau 0 equations was developed from 
combined June 1976 and June 1977 data and verified with the 


dependent data. 


Coefficient Predictor 
-303.043 
+ 1.165 EHF 
+ 35 BVISR 
e اانا‎ OB PS 
+ 4.627 SEHE 
SI ASTDR 
- 2.609 VCMP 1 
ee 
Dependent Verification: Percent Correct = 39 
Skill Score = „13 
Category 1 T LIF IV V Vi VII VELER TX X 
Bias O CO0 OL 201205 —09 .4 I mo 00 


This equation shows some improvement over the l0CATA 
equation in R^ value, however the percent correct is slightly 


lower and the Heidke skill score is the same. 
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es Caregory scheme:  >CAT 


Deriving a regression equation with fewer categories 


should yield better results due to partial compensation of 


observer error. 


In this case, 


five categories are used 


which correspond to the probabilistic five-category scheme 


of Aldinger (1979). 


Category Observed 
Number Visibility Codes 
1 9079 1792 

51 93794 

Ji 337736 

IV 57 

V 90799 


Visibility Value of 
Range (km) Predictand (km) 
OS 505 
OSTEO «72.0 25 
Zero <10.0 6:0 
lO.O Eo <20.0 15.0 
>20.0 SERO 


The Tau 0 equation was developed from combined June 


1976 and June 1977 data, and verified using both the dependent 


June data and independent data from July 1979. 


Coefficient 
ZO 
797.035 
02? 
e 7 
174,280 
+ 0,944 
2000 215 


Predictor 


EHE 
BVISR 
PS 
SERF 1 
ASTDR 


VCOMP 
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Dependent Verification: Percent Correct 


Es Score = .17 
Category I TI PIT IV V 
Bias .02 92 .47 2 2 17.05 
Independent Verification: Percent Correct = 42 
SUME re = .17 
Category 1 ET EIT IV V 
Bias -03 -02 :25 297 .49 


It is to be noted that the variables selected are the 
same as those selected in the two ten-category schemes with 
the exception that in this scheme VCOMP was selected instead 
of VCMPl. The 5CAT scheme shows an increase in skill score 
as expected, and the percent correct also increased. Bias 
values here are not much better than those for LOCATA and 
lOCATB except for category V of the dependent verification 
and category IV of the independent verification, both of which 


show values approaching unity. 


B. REGRESSION EQUATIONS 

The ultimate goal is to forecast, not just analyze, visi- 
bility. Therefore, using the July 1979 data set and a new 
set of parameters which included prognostic predictors, new 
equations were developed using the 5CAT scheme. First a new 
equation for Tau 0 was derived, then forecast-interval equa- 


tions for Tau 24 and Tau 48 were developed. The parameter set 
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used for these equations is given in Appendix A, part 2, 
All three of the following equations were verified using the 
dependent data and also verified independently with data from 


August 1979. 


MO0=n! Diagnostic Equation:  5P0O0 


Coefficient Predictor 
US 
EEO G7 EHF 00 
+ 0.488 BVISR 
- 98 FTER 00 
+ 3.048 SEHE IS 2 
Rf = .30 


The two-digit number after some of the predictor 
parameters indicates the time interval from which the 
parameter is derived. Those predictors without such a number 


are available at the analysis time only. 


Dependent Verification: Percent Correct = 42 

SIS core = .18 

Category 1 LI III IV V 

Bias 202 02 . 90 2:027 1207 

Independent Verification: Percent Correct = 5l 

51-1 11500126 = .21 

Category I II IIL TY V 

Bias 207 802 299 22:00 1-10 
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The m value and verification of equation 5P00 is 
better than the verification of the 5CAT equation due to the 
consideration of more parameters in the July 1979 data set 
Eun the combined June 1976 and June 1979 data sets. The 
bias values are not much different, except for category III 
which shows improvement. It may be noted that all selected 
parameters but one are from the analysis time which seems 
consistent with the nature of the Tau 0 equation. 

An interesting fact is that the independent verifica- 
tion of 5P00 yields better values than the dependent verifica- 
tion. This is, in part, due to the fact that the independent 
data contains a higher percentage of observations in those 
high visibility categories which the equation estimates best. 
In addition the dependent data comes from a large enough 
sample of synoptic conditions that the regression equation 
could score higher when applied to independent data, which 
by chance includes a larger number of those synoptic situations 
best handled by the equation. 


P 24—hr Prognostic Equation:  5P24 


Coefficient Predictor 
SS 
E T07 EHF 24 
+ 0.440 BVISR 
E O 002 RHRX 
- 7.418 FTER 24 
R^ - .30 


34 





Dependent verification: Percent Correct = 42 


Skill Score = ‚16 
FZregory I TI TILI EN V 
Bias 10 . 08 56 2.26 16 
mdiesendentsverifieation: Percent Correct = .52 
Skill Score = .20 
Sage gory I TI ITI IV V 
Bias .04 UN “CT 10,92 2817 


There is a deterioration in Rd value when 5P24 is 
compared to 5P00, as one might expect. The percent correct 
is similar for both equations, but the Heidke skill score for 
5P24 is slightly less than for 5P00. Here again, as in 5P00, 
the independent verification is better than the dependent 
verification. 

It is to be noted that variables from Tau 24 have 
entered the 5P24 equation, which is consistent with the 
nature of a Tau 24 equation. 


3. 48-hr Prognostic Equation:  4P48 


Coefficient Predictor 
- 4.160 
10.5390 EHF 36 
1280/5955 BVISR 
=12. 631 FTER 48 
70.033 EHF 00 
+ 0.003 RHRSO 
n 0.160 MBVIS 48 
4? 5 
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Dependent Verification: Percent Correct = 42 
Skill Score = „13 
Category nh TI TIE IV V 
Bias TOT SURE 229 2203 1.40 
Independent Verification: Percent Correct = 52 
Skill Score 6 
Category d II TIT IV V 
Bias . 00 wou “20 1972 132 


Here the ES value has deteriorated somewhat from the 
5P00 and 5P24 cases. The percent correct is the same for 
equations at all three time periods, but the Heidke skill 
score in 5P48 is worse than that for 5P24 and 5P00. Overall 
the bias values for 5P48 are worse than for both 5P00 and 5P24. 
Once again the independent verification is better than the 
dependent verification. 

It is to be noted that two Tau 48 hr predictors have 
entered the equation. However, there is also one TAu 36 hr 
predictor and three Tau 00 hr predictors. The predictor 
BVISR shows up in 5P48 as well asin 5P00 and 5P24. BVISR, 
which itself is a parameterized visibility, can be considered 
an indicator of the persistence of marine visbility regimes 


through 48 hours. 


C, PROBABILISTIC VS. CATEGORICAL APPROACH 
Aldinger (1979) used the 5CAT scheme outline previously 


and developed regression equations for the probability of 


36 





ereurrence of each category. Then, using the notion of 
threshold probability, the most-likely category was determined. 
For comparison, an equation was developed by the categorical 
method of this study considering only those predictor parameters 
used by Aldinger. All equations were derived from the com- 
bined June 1976 and June 1977 data and were verified dependently. 


I Probabilistic Equations [Aldinger, 1979] 


Category Equation 
I VISPROBT T= 36067262 > 1.647 SEHE + „2839 RHR 


= „369 PS + „401 VCOMP 


B^ 2 .13 
II VISPROB - 738.837 - .264 EHF - .746 PS 
+ „555 RHR - 1.689 SEHF 
in. cs A 
III VISPROB = 266.075 + .303 VVWW - .256 PS 
+ .247 RHR + .313 RHX 
R^ = .5 
IV VISPROB = -278.669 + .365 SEHF - .643 VCOMP 
+ .431 VVWW + .333 PS 
R* = .9 
V VISPROB = -693.510 + 3.633 EHF + .767 PS 
- .709 VCOMP - .352 RHR 
2ج‎ - .1 
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VISPROB is the probability of occurrence of the category 


for which the equation is derived. 


Dependent Verification: Percent Correct = 32 
Skill Score = .13 
Category T geb Ju IV V 
Bias .04 1.58 11,0 203 0.40 


2. Categorical Equation 


Only one categorical equation was derived whose 
visibility value (VIS) determines the visibility category 
by selecting that category to which VIS belongs. 

IS 502,35 + 2175 BHF + .339 PS - .254 RHR 

EDO SENE 


RZ = .24 


Dependent Verification: Percent Correct = 43 
Skill Score = „14 
55 900217 1 TI TTI IV V 


Bias 202 Qu 28 2-068 T3 


Comparing the two approaches shows that the cate- 
gorical approach yields a higher percent correct and a 
slightly higher skill score. However, except for category 
V, the biases are worse for the categorical scheme As might 
be expected both methods use similar predictor parameters. 
SEHF, RHR, PS and EHF are common to both. 

3. NPS Awards Matrix 
Aldinger (1979) developed an awards matrix which 


when applied to the verification matrix (Appendix E) of a 
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5-category scheme gives some credit to near successes. The 
Techniques Development Laboratory (TDL) of the National 
Weather Service has also used an awards matrix, but of a 
different nature, which does not give full credit to all 
correct visibility estimates [National Weather Service, 1973]. 
The NPS awards matrix does give full credit to all correct 
estimates. All quantities of a verification matrix are 
multiplied by the corresponding percentages in the awards 


matrix shown below. 


OBSERVED Estimated Category 
CATEGORY 13 EI TIT IV M 
I 100 80 0 0 0 
mer 80 100 25 0 0 
ET 0 25 100 25 0 
IV 0 0 25 100 13 
V 0 0 0 15 100 


The verification results, after applying the awards matrix, 


are as follows: 


e Drele Approach: Percent Correct = 60 


Skill Score = .27 
Categorical Approach: Percent Correct = 63 
Skill Score - .2 


In both cases percent correct increases markedly. 


However, for the probabilistic approach the skill score doubles, 
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while for the categorical approach the skill score decreases. 
This shows that the probabilistic approach forecasts near 
successes much better than the categorical approach, thus 


enhancing its usefulness. 


D. PREDICTAND TRANSFORMATIONS 

Generally the relationship between an atmospheric pre- 
dictand and the predictors is not linear. This can lead to 
less than desirable results when multiple linear regression 
as used. Non-linear regression may be used to overcome this 
problem, but the increased computational time involved usually 
precludes its use. Another method used to solve the non- 
linear problem is to transform the predictand to a form which 
then relates in a more linear manner to the predictors. 

Using a limited number of parameters several transforms 
were tested on the LOCATA scheme, using July 1976 and July 
1977 data. The relative values of R^ produced using each 


transform are shown below. 


Predictand e 

VISIBILITY (VIS) .230 
Log, y (VIS) 243 
1/VIS .037 
(1/vıs)? 2011 
vıs!/2 570 
vis? re 
visi/4 E 
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It can be seen that the Be value for several of the 
transformed predictands was higher than the ae value for 
Siemon -transftormed visibility predictand, though the 
increase was not large. 

However, the real test is how well an equation with a 
transformed predictand verifies. So the equation derived 
with the cube root of visibility as the predictand, which 
yielded the highest m value, was scored against the equation 
with the non-transformed predictand. 

Beedictand = visibility. 


S9 


Dependent Verification: Percent Correct 


Skill Score = .14 


Predictand = visibility 


DE me na er | Cartons Percent Correct = 27 
Skill Score = -,0l 


The results show that the transformed predictand yielded 
worse scores than the unmodified visibility predictand. 

This is a surprising result in view of the relative R value: 
It may, in part, be explained by the fact that there was an 
uneven distribution of visibility observations between cate- 
gories, with a heavy weighting toward higher visibility cate- 
gories. Time limitations, however, did not permit examining 
this further, and all other research was conducted using the 


non-transformed predictand. 
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E. WEIGHTED LEAST SQUARES 

In this study the data distribution is such that most 
observations occurred in the higher categories, in particu- 
lar category 98. The result of this is a regression equation 
that fits the higher visibility categories better than the 
lower visibility categories. As a result, low visibilities 
are poorly estimated. 

The technique of weighted least squares was applied in 
an attempt to alleviate this problem. The goal was to weight 
more heavily the lower category cases in relation to those 
in the higher categories so that the resultant equation would 
increase skill in estimating poor visibilities. 

The BMDP programs [UCLA, 1979] allow case weights to be 


applied. The weighted least squares technique minimizes 


2 2 
IR Sr ET 


where, 


E is the case weight for case j 


ve is the observed visibility for case j 


^ 


Y is the regression estimate for case j. 


Normally the weight for each case should be inversely 
proportional to the variance [Daniel, 1971], but any number 
of weighting techniques may be tried. In this study two 
sets of case weights were tried and applied to the schem of 


lOCATA. 
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The first scheme (WLS1) weighted each case with a weight 


equal to the inverse of the predictand value, as follows. 


For cases of The predictand And the case 
observed code value (km) is weight Ml is 
90 025 17.2025 
31 25 197.125 
92 35 14.35 
93 75 1% 21.3 
94 O 1.5 
95 t) Ip». 0 
96 7900 177-0 
97 lS 0 W150 
98 35.0 1032.09 
99 MORO 75.0 


The resultant equation derived from combined June 1976 
and June 1977 data (not given here) was verified dependently 


with the following results. 


Percent Correct = 7 


=e OL 


Skill Score 


Obviously, this is a poor weighting system. The ES value 
is very low and the scores are predictably poor. 

For the second scheme (WLS2) a more reasonable set of 
weights was used. The variance was computed for each cate- 


gory from the unweighted equation of 10CATA. Then the weight 
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for each case in a particular observed category was set to 


the inverse of the square root of the variance of the observed 


category. 
For Cases of The Predictand And the case 
Observed Code value (km) is weight (wa) is 
90 025 0052 
91 25 . 0603 
92 35 . 0661 
93 275 -0615 
94 5 20:702 
95 3.0 .0700 
96 7.0 .0754 
57 15:0 . 0941 
98 35.0 80925 
99 pe» «0 .0242 


(Each code group corresponds to a category in the 10CATA 


scheme.) 


The case weights shown here are somewhat contrary to what 
might be expected. It would seem that the variances of the 
higher categories would be larger than those of the smaller 
categories, if for no other reason than the fact that the 
visibility ranges of the higher categories are greater. If 
this were true the case weights for the higher categories 
would be smaller than for the lower categories. However, 


the weights shown here generally increase with an increase in 
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category, with the exception of category X (code 99). This 
result is due to the fact that the regression equation esti- 
mates those categories best which contain the highest number 
of observations, namely the categories containing codes 97 
and 98. 

A comparison of dependent verification between the equa- 


tions of 10CATA and WLS2 shows very little difference. 


Scheme R? Percent Correct Skill Score 
BEOCATA 125 40 3 
WLS2 23 40 2 

2 


EH DEFLATION OF R 

According to theory, if a regression equation perfectly 
fits the data from which it was developed the explained 
variance, should equal 1.0. However, it appears that 
due to the nature of the categorical schemes in this study 
a limit was placed on the maximum R* that it was possible to 
achieve. This particular limit is related to the fact that 
each predictand value was assumed to be the midpoint value 
of the observed category, thus providing discrete visibility 
values. However, the regression equation gives continuous 
visibility values which are then used with the assigned pre- 


dictand values to determine R^. 


In one experiment, to demonstrate the deflation of R^, 
a regression equation of the form of l0CATA scheme was 


developed. Then using the dependent data, the equation was 
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used to compute visibility values, Vi. 


1 cT II 


Symbolically: V; = Aj + B,X,; je 2 


1 il 


where, 


V visibility 


x's independent predictors. 


These Vi values were used as substitutes for the original 
visibility observations. Next, using these Vi values, a new 
predictand, Vi', was derived by re-setting the Vi value to 
the midpoint of the category to which Vi belonged, giving 
Vite Finally, a second regression equation was developed 


using the Vif as predictand values to yield an equation of 


the form 


Tf = 
E A, + BJX]; + CoXoj 0 RA 


It can be seen that if the continuous values, Ver had been 
used as the predictand the second regression equation would 
be identical to the first one and have an R? value of 1.0. 
However, because the predictand, Vir used to develop the 
second equation has discrete values as defined by the cate- 
gorical scheme, the second equation is not identical to the 
first; and the R^ value is approximately 0.7, using V, ' as 
the observed values. 

It is believed that the x value of 0.7 rather than 1.0 


is the maximum value achievable in the LOCATA scheme with a 
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perfect equation, due to the method of defining the pre- 
dictand used in this study. The other categorical schemes, 
of course, have a similar R? mares 

Ene drop of Be ال ا‎ UO TG U: can De demonstrated 
by schematic graphs. Assuming that the observed visibility 
can be expressed perfectly by a regression equation, for 
which R? = 1.0, then the graph below is the result. As 


the continuous regression-estimated visibility increases 


the observed visibility increases continuously also. 
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Regression Equation 
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However, the observed visibility is not given as a 
continuous variable. Rather the visibility observations 
are given as ranges or categories, and the visibility 
predictand is defined as the midpoint of the observed 


range, which is demonstrated schematically below. 
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Visibility from Regression Equation 


The schematic above shows a step function relationship 
which indicates that as the continuous regression- 
estimated visibility increases within each categorical 
visibility range (given by roman numerals) the observed 
visibility remains constant. 

The regression-estimated visibility values have not 


changed from the first schematic to the second but the 
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verifying "observed" values have changed from continuous 
to discrete values. All observed values below a categorical 
midpoint value have been increased, and values lying above 
a midpoint value have been decreased. 

The deterioration of R* which results from the second 
case can be seen by noting the deviation of values along 
the discrete observed visibility step function from the 


continuous observed visibility line as shown below. 


| deviation 
(discrete -continuous) 


Observed Visibility 


c 


Visibility from Regression Equation 


In another experiment, an attempt was made to compute 
the pé value for the LOCATA equation without the hindrance 
of the problem just described. The BMDP programs compute 


R using the continuous regression-produced visibility 
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values and the discrete observed values. A separate program 
was developed to compute R* by first re-setting the continuous 
regression values of 10 CATA to the midpoint values of the 
categories to which they Pe Then, using the discrete 
predictand values, a new R* was computed. In this case dis- 
crete values are used for both the observations and the 
regression estimates. The BG value computed in this way is 
.3l as compared to .25 computed by the BMDP programs. All 

Ee values previously shown in this study were computed by 

the method used in the BMDP programs. 

The maximum R* value of approximately 0.7 as found by 
experiment for the LOCATA scheme may be compared to the R* 
value of .31 which the 10CATA equation yielded. The differ- 
ence between the two uc values of approximately 40$ can now 
be attributed to errors in the observations and numerical 
MOP's and the non-linear relationship between visibility and 


associated meteorological parameters. 


G. DISTRIBUTION PROBLEM 

The distribution of observations among synoptic codes for 
the combined June 1976 and June 1977 data set is shown below. 
It can be noted that the highest three categories contain 
66% of the observations, and the highest four categories 
contain 79% of the observations. The observation distribu- 
tions are similar for the July 1979 and August 1979 data 


sets. 
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Code Number of Percent of 


Group Observations Total observations 
90 75 0.8 
51 2:39 206 
22 400 4.4 
33 740 Su 
94 166 1.8 
95 327 3.6 
96 1125 1.208 
p Dorn 210 
93 3642 DOM 
99 495 5.4 


This fact tended to tune all the regression equations to 
the high categories, such that high categories were estimated 
relatively well by the regression equations and low visi- 
bility categories were estimated poorly. This is somewhat 
contrary to what is desired, since forecasts of low visibility 
are very important operationally. 

The probabilistic approach does not have a similar dis- 
tribution problem, since one regression equation is developed 
for each visibility category and depends only on the observa- 


tions of a single category. 


E BETA VISIBILITY 
The beta visibility was previously described. Its compu- 


tation is given in Appendix B.3. Beta visibility is not only 
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a parameter for use in visibility regression equations but 
itself yields a value of visibility which may be of use. 
This section attempts to quantify its usefulness. 

The BMDP programs were used to compute a correlation 
coefficient between the predictand and the various forms of 
the beta visibility parameter. It is to be noted that the 
visibility predictand is not a directly observed visibility 
value, but rather it is the midpoint value of an observed 
visibility range. The correlation coefficients, R, between 
the various forms of the beta visibility parameter and the 
visibility predictand of the 5CAT scheme are given in the 
following table. A comparison of maximum, minimum and mean 
values is also given. These statistics were derived using 


ner July 1979 data set. 


Comparative Statistics and Correlation to the Visibility 
Predictand (VIS) at Tau 0 hr 


Maximum (km) Minimum (km) Mean (km) E 
VIS (Tau 0) 33.0 025 192 1-00 
BVISR 46.9 056 14.3 0.43 
BVISX 5 O79 OS 0:09 


Comparative Statistics and Correlation to the Visibility 
Predictand (VIS) at Tau 0+24 hr 


Maximum (km) Minimum (km) Mean (km) R 
VIS (Tau 24) 350 0525 190 00 
BVISR 48.7 Oreo! 14.3 Oe Suk 
BVISX 259 DES 20160 O10 
MBVIS 24 44.4 1.68 1752 0205 
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Comparative Statistics and Correlation to the Visibility 
Predictand (VIS) at Tau 0+48 hr 


Maximum (km) Minimum (km) Mean (km) By: 
VIS (Tau 48) 35:0 0.25 oo O0 
BYTSR S21 0.42 14.3 0.24 
BVISX Su 062 20.0 005 
MBVIS 48 30. l 2.14 15.4 0.02 


It should be noted that in the table the analysis- 
time parameters BVISR and BVISX are compared to the 
predictand at all three time periods. The table shows 
that the maximum, minimum and mean values of all the beta 
visibility parameters are similar to the corresponding 
values of the visibility predictand at each time period. 
BVISR shows a higher correlation to the predictand than 
BVISX at all time periods, though the correlation of both 
parameters to the predictand worsens with time. Both the 
analysis-time parameters BVISR and BVISX show higher 
correlation to the predictand at Tau 24 hr than the 
prognostic-time parameter MBVIS 24. The same is true at 
Tau 48 hr when comparing BVISR and BVISX to MBVIS 48. 

The following clarifies the reason for the slight 
differences in maximum, minimum and mean values for the 
same parameter at different time periods. The Tau 24 hr 
data includes values from the first day of August (i.e. 


up to 24 hrs after the last day of the July data set), and 
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omits values from the first day in July. In like manner, 
the Tau 48 hr data includes the first two days of August 
and omits the first two days of July. Thus the data set for 
each time period is slightly different. 

In addition, a skill score was computed for BVISR and 
BVISX by determining the code group to which the computed 
beta visibility belonged, and comparing that to the observed 


code groups in the combined June 1976 and June 1977 data. 


Heidke Skill Score Percent Correct 
BVISR 0.10 33 
BVISX 0.07 31 


It can be concluded by these results that although beta 
visibility is a useful predictor parameter for regression 
analysis, it has quite limited skill when used to estimate 


visibility by itself. 


I. COMMENTS ON EXPLAINED VARIANCE 

The total explained variance, me, of a multiple linear 
regression equation is a measure of how well the dependent 
variable (predictand) can be approximated by a linear com- 
bination of independent variables (predictors). The higher 
the value of R^, the better the approximation is. A perfect 
linear relationship results in an R? value of 1.0. However, 


it should be noted that R* indicates only how well a given 


equation will estimate a given predictand if one uses the 
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method of least squares. This method results in a regression 
equation which minimizes the value of the sum of squares 

of the estimate errors (estimate error = estimated value minus 
observed value). An equation with a given R* will not 
necessarily provide a better estimate of the predictand than 
an equation with a smaller R* when evaluated by some method 
other than least squares. An entirely different situation 
may occur if one applies the derived regression equation to 
independent data. Though the original equation may be a 

good fitting equation for the dependent data (by the least 
squares criterion) it may be a poor fit for the independent 
data, especially if the number of cases is small. In this 
study the sample size of over 4000 cases is large enough that 
a drastic drop in estimation ability is not to be expected 
when independent data are applied, however some deterioration 
Ew encountered. 

Also, as additional predictors are entered into an equa- 
tion by the stepwise process the R* value will increase, but 
an equation with fewer predictors and a lower n? May, in fact, 
provide a better estimate when applied to independent data. 
This is so, since as more variables enter into an equation, 
it becomes more likely that the equation will reflect relation- 
ships unique to the dependent data. Thus extra variables 
may degrade an equation when scored on independent data [Air 
Weather Service, 1977]. Of course, the application of inde- 


pendent data may also show an improvement in scores due to 
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the peculiarities of a particular data set. However, some 
form of truncation method should be used to limit the number 
of variables in an equation such as was done in this study. 
An experiment to demonstrate the relationship of score 
to number of predictors in the equation was performed, using 
the regression results of the 5CAT scheme. Truncating the 


SCAT scheme at different steps yielded the following. 


= Dependent Data Independent Data 


Skill Score SECOrrece Skill Score % Correct 























40 21223 Dr 





.4 





208 








42.7 z AS 41. 





. 9 











SS 44.0 20579 42. 














43 . 





Al 43.2 73 











167 43 .8 179 42. 








41. 






165 





44 





. 4 20 











. 166 44.4 #156 41. 













.174 186 > 40.9 





It can be seen that after a certain point the direct 


relationship between I and skill becomes obscure. In 
this study the equation for the 5CAT scheme as described 
in the text was truncated after the sixth step, for at 
the seventh step the R failed to increase by a rounded 


value of 13. 
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It is encouraging to note that the results above show 
that percent correct and skill score do not substantially 
decrease when independent data is applied compared to when 
dependent data is applied. In fact, the skill score is 
relatively better in the former instance for the first 


five steps. 


DISCUSSION OF ERRORS 

It is believed by the author that the techniques used 
in this study would yield equations of high operational 
usefulness if it were not for various unavoidable errors. 
Linear regression assumes, for example, that all predictand 
values used are errorless. This is far from true here. 
Observer error in estimating visibility at sea is relatively 
high, due mostly to a dearth of visibility markers at sea 
and also due to the fact that many ships transmitting 
synoptic reports may have observers with little or no 
observational training and/or experience. 

Errors also enter into the Model Output Parameters, 
which are only as good as the numerical models from which 
they are generated, analyses being better than prognosis. 
The method used to interpolate the MOP's to the synoptic 


ship positions also adds error to the scheme. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The categorical approach used in this study yielded 
Visibility equations which have comparable skill both at 
analysis and prognostic times which iS a promising result. 
However, the actual skill of the equations is relatively 
poor and not operationally useful at this time. The 
reason for this is believed to lie inherent in the errors 
of visibility observations, the non-linear relationship 
between the predictand and the predictors, and the 
numerically generated MOP's. The future promises much 
improvement due to new statistical techniques, improved 
numerical models and the identification of more air/ 
ocean parameters with a known relation to visibility. 

The comparison of the probabilistic to the categorical 
approach indicates that the probabilistic approach holds 
more promise, at least partly due to the fact that the 
categorical approach is hindered by the uneven distribution 
of observations. The probabilistic approach seems to 
estimate near successes better than the categorical 
approach. 

Parameters found to be most highly related to visibility 
in the regression equations are: evaporative heat flux, 
beta visibility, sea level pressure, sensible plus 


evaporative heat flux, air/sea temperature difference, 
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meridional component of the wind, relative humidity 
parameters and FNOC's fog probability parameter. 

The following recommendations are offered for future 
research: 

l. Test new parameters in relation to visibility, 
such as some type of visibility persistence parameter, 
more interactive, modified and binary parameters, and a 
climatological parameter now being developed for the 
North Pacific by the National Climatic Center. 

2. Investigate further the techniques of weighted 
least Squares and transformation of the predictand to 
relate more closely to the non-linear nature of the 
problem. 

3. Stratify the data with respect to critical values 
of geography and to various MOP's. 

4. Investigate the use of discriminant analysis to 
estimate visibility. 

5. Stress the probabilistic approach over the 
categorical approach, and in particular, expand the 
work of Aldinger [1979] to include additional parameters 


and prognostic equations. 
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Part l. 


APPENDIX A 


PREDICTOR PARAMETER DESCRIPTIONS 


This part consists of all predictor parameters 


considered for use in the analysis-time equations 


developed from the combined June 1976 and June 1977 data 


set. 


NOTES : 


a 


[*] 


=] 


SYMBOL 


Denotes those predictor parameters that 
repeatedly were selected early by the stepwise 
regression thereby implying their relatively 
strong relationship with visibility. 


Denotes those predictor parameters that only 
occasionally or never were selected early by 
the stepwise regression, but may be useful in 
future studies. 


Denotes those predictor paramters that seemed 


comia werde OF nO relation to visibility in 
this study. 


DESCRIPTIVE NAME UNITS 


A. Analysis Parameters (FNOC Mass Structure Model) 


PS 


TAIR 


EAIR 


T925 


TSEA 


Sea-level Pressure [**] (mb) 
Surface Air Temperature [*] (2C) 
Surface Vapor Pressure [*] (mb) 
925 mb Air Temperature [*] (€) 
Sea-Surface Temperature [*] (°C) 
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B. Prognostic Parameters 


TX Surface Air Temperature [*] 
Derived from surface air and potential 
temperatures, boundary layer depth, 
upper-level winds extrapolated to 
surface, air density, drag coefficient, 
gustiness factor 


EX Surface Vapor Pressure [*] 
Derived from model's mixing ratio 


SOLARAD Solar Radiation [*] 
Calculated absorption of incoming 
short-wave (solar) radiation. 
(postive downward) 

EHF Evaporative Heat Flux [**] 
Derived using air density, drag 
coefficient extrapolated winds, 
and mixing ratios. 

SHF Sensible Heat Flux [*] 
Recovered from SHF = SEHF-EHF. 
Originally derived by FNOC using 
drag coefficient, extrapolated winds, 
surface air temperature, TX, 
density and constants. 

SEHE Sensible Plus Evaporative Heat Flux [**] 
SERF = SHE+EHEF 

EHF Total Heat Flux [*] 


THF = SEHF-SOLARAD+LW, 
where LW is the heating due to long- 
wave (terrestrial) radiation. 


C. Marine Wind Model (FNOC) 


VVWW Marine Wind Speed [*] 

Marine Wind Direction 

This variable was not used as a 
predictor parameter, but rather 
to derive other parameters. 


(DDWW) 


D. Derived Parameters 


UCOMP Zonal Wind Component [*] 


UCOMP = -VVWW sin(DDWW-10) 
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and empirical constants. 


(FNOC Primitive Equation Model) 


Oo 


(mb) 


(sea 
cm? /hr) 


(gear, 
cm2/hr) 


(gcal/ 
cm2/hr) 


(gcal/ 
cm2/hr) 


(gel! 
cm2/hr) 


(kt) 


(deg/10) 


(m/sec) 





VCOMP 


CAPU 


CAPV 


THETAX 


THETAR 


STABX 


STABR 


ASTDX 


ASTDR 


ADTSEA 


ADTX 


ADTAIR 


AASTDX 


AASTDR 


Meridional Wind Component [**] 
VCOMP = -VVWW cos (DDWW:10) 


I Directional Wind Component [*] 


CAPU = -UCOMP - sin(LNGA) 

-VCOMP ° cos (LNGA) 
[Haltiner, 1971], where 
ENGA = -10 - (I,J point longitude). 


J Directional Wind Component [*] 
CAPV = VCOMP ° cos (LNGA) 
-VCOMP - sin (LNGA) 
[Haltiner, 1971], where 
CNGA = =]0= (I,J point longitude). 


Potential Temperature X [-] 
Derived using PS, TX. 


Potential Temperature R [-] 
Derived using PS, TAIR. 


Stability —| 
Derived using [THETAX - 
(THETA from T925)]/(PS-925]. 


Stability R [-] 
Derived using [THETAR - 
(THETA from T925)]/(PS-925). 


Air-Sea Temperature Difference X [**] 
ASTDX = TX-TSEA 


Air-Sea Temperature Difference R [**] 
ASTDR = TAIR-TSEA. 


Advection of TSEA [*] 
See Appendix B.l. 


Advection of TX [*] 
See Appendix B.l. 


Advection of TAIR [-] 
See Appendix B.l. 


Advection of ASTDX [-] 
See Appendix B.l. 


Advection of ASTDR [*] 
See Appendix B.l. 
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(m/sec) 


(m/sec) 


(m/sec) 


(°K/mb) 


(°K/mb) 


(“C7 EE) 


(CRE) 


(CERT) 


) E 


(SC ANE) 






Veep iape 






۰ nm 


RHR Relative Humidity R [**] ES 
See Appendix B.2. 

RHX Relative Humidity X [**] (%) 
See Appendix B.2. 

E. Interactive and Modified Parameters 

RHRX = RHR + RHX [**] 

RVCOMP = RHR - VCOMP [-] 

RHRPS = RHR + PS [-] 

RASTDX = RHR + ASTDX [**] 

RSEHF = RHR ° SEHF [-] 

PDSQ = (PS-1014.8)* [-] 

PSRHX = PS - RHX [-] 

PSSEHF = PS - SEHF [-] 

PASTDX = PS + ASTDX [*] 

PSVCMP = PS + VCOMP [-] 

VSEHF = VCOMP + SEHF [-] 

EHFADT = EHF + ADTAIR [-] 

ESEHF = EHF + SEHF 

EXEAIR = EX * EAIR [-] 

SEVCMP = SEHF - VCOMP [-] 

SEADTX = SEHF + ASTDX [-] 

SERHX = SEHF * RHX [-] 

ASTDRX = ASTDR + ASTDX [*] 

UVCOMP = UCOMP » VCOMP [*] 

CAPUV = CAPU - CAPV [*] 

TARSEA = TAIR * TSEA [-] 

TXAIR A A |] 
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SEHFSO = SEHF ° SEHF [-] 


PRESO = EHF - EHF [-] 

RHRSQ = RHR + RHR [(**] 

RHXSQ = RHX + RHX [*] 

VCMPSO = VCOMP - VCOMP [-] 
CAPUSQ = CAPU - CAPU [*] 

TSEASQ = TSEA - TSEA |-] 

ASDXSO = ASTDX ٠ 2511015 [^*] 
ASDRSO = ASTDR ° ASTDR [*] 
ADSESQ = ADTSEA - ADTSEA [-] 
PSSQ = PS ° PS [-] 

SREHF SGUaré FOOL Of EHF [7] 
SRPS SaUars FOOL ot PS 1%] 
SRASTR Square root of ASTDR [-] 
SRASTX Square oot of ASTDRX [=| 
SRSEHF Sanare I OOt Of SEHR [*] 
SRRHR Seuare OOC Of RHR [=] 
SRRHX Scuüare 2500 of RHX [2| 
SRCAPU 56002226 2068 of CAPU 1=] 
SRTSEA Square root of TSEA [-] 
SRVCMP Square root of VCOMP [-] 
SRASEA Square root of ADTSEA [*] 


F. Binary Parameters 


EHF1 if EHF «1.75 or EHF » 8.75; EHFl = 0.0 [-] 
if 1.75 « EHF « 8.75; EHFl - 1.0 
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EHF2 if EHF < 3.34; EHF2 = 0.0 [s] 
if EHF > 3.34; EHF2 = 1.0 
EHF3 if EHF < 0.0; EHF3 = 0.0 [-] 
ا‎ = 0.0,  CEHES = 1.0 
PS1 if PS < 1000 or PS > 1030; PS1 = 0.0 [-] 
jo 28 100, PS1 = 1.0 
PS2 11 25 > 1014.3: PS2 = 0.0 [-] 
E ES IOLA SG ES2 = 1.0 
RHR1 if RHR < 60; RHR1 = 0.0 [-] 
Bere CO FIRI = 1.0 
RHR2 if RHR < 83; RHR2 = 0.0 [-] 
E EIR. 6: RHR2 3 1.0 
SEHF1 if SEHF « 0.0; SEHFl = 0.0 [**] 
MMC OO; SEHEL = 1.0 
ASDX1 FASTO 30 0 -SASDX] - 0.0 [=] 
if ASTDX > 0.0; ASDX1 = 1.0 
ASDR1 if ASTDR < 0.0; ASDRI s 0.0  [-] 
TE ASTON 0<0; ASDRI = 1.0 
VCMP 1 i00 200: VEMPL = 0.0 £x[**] 
CONE 0.0; VeMPl = 1.0 
UCMP1 E UCONN 0.0, UCMPL = 0.0 [=] 
jM COMPE 0-0; UCMPL = 1.0 
STABX1 if STABX < 0.0; STABX1 = 0.0 [-] 
ESTADO 0) STABKI = 1.0 
STABRI if STABR « 0.0; STABRl = 0.0 [-] 
ASTAR - 0.0; STABRI = 1.0 
G. Other Parameters 
ETER FNOC Fog Probability Parameter [**] (3) 
BVISR Beta Visibility Parameter R [**] (km) 
See Appendix B,3. 
BVTSX Beta Visibility Parameter X [*] (km) 


See Appendix B,3. 
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Part 2. This part consists of all predictor parameters 
considered for use in the analysis-time and forecast- 
interval equations developed from the July 1979 data. 
In this list some parameters not found useful in the 
June regression runs were eliminated, but additional 
parameters which were available for the July data set 


were added. 


A. Predictors used to develop equations both from June 
and from July data (described in Part 1) 


(1) Parameters available for Tau 00, 12, 24, 36 


and 48 hr 

PS T925 TX 

EX EHF SHE 
SEHE 1HE VVWW 
UCOMP VCOMP RHX 
EHF2 SERE IL VCMP 1 
FTER UVCOMP 


(2) Parameters available for Tau 00 hr only 


TAIR EAIR TSEA 
ASTDX ASTDR RHR 
ASTDRX ASDXSQ RASTDX 
RHRX RHRSQ BVISR 
BVISX 
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B. Additional variables available in the July 1979 


data set 

SYMBOL DESCRIPTIVE NAME UNITS 

CLIMO National Climatic Center (3/100) 
Fog Frequency Climatology [*] 

SSANOM Sea Surface Temperature Anomaly [*] me) 
Available at Tau 00 hr 

W925 oenina Component at 925 mb [*] (KE) 
Available at Tau 00, 12, 24, 36, 48 hr 

2925 V Wind component at 925 mb [*] (kt) 
Available at Tau 00, 12, 24, 36, 48 hr 

E925 Vapor pressure at 925 mb [*] (mb) 
Available at Tau 12,24,36,48 hr 

GGTHTA Front Location Parameter [*] Uo > 
Available at Tau 00, 12, 24, 36, 48 hr (100 km) ) 

NCLOUD Total Cloud Cover 1*] (tenths) 


Avallable at Tau 00, l2, 24, 36, 48 hr 


MBVIS Modified beta visibility [**] (km) 
See Appendix B.3 
Available rat Tau 12, 24, 36, 48 hr 


RASTDR = RHR * ASTDR [*] (°C $) 
Available at Tau 00 hr 


H510 0006ل‎ 2 AA (cm) 


D-value thickness 
Available at Tau 00, 12, 24, 36, 48 hr 
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APPENDIX B 


MISCELLANEOUS PARAMETER FORMULATIONS 


l. Advection Parameters 

All advection parameters use the following general 
formulation. 

For the advection of a quantity (Q) the formula 


ADQ = ES * VQ, was used in the finite difference form: 


RMAP 


DM ], 


ADQ = - 





comer GY O — 97.3) 


J JE 


where RMAP = (1 + sin 60)/(l + sin (latitude) ) 


LS ano + sini 60) 1/31. 205 


and DM 
(31.205 = number of grid mesh lengths, pole to equator, 


n Elle FNOC I,J grid). 


2. Relative Humidity Parameters 
The thermodynamic equation for calculation of 
Saturation vapor pressure, known as the Clausius-Clapeyron 


equation is given as 


E L(T)/RT“ : (m 


where 
R = specific gas constant for water vapor 


(0.461 joule gt ek 1) 
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T = temperature (°K) 


L(T) = latent heat of vaporization of water 
(joule g ") 
ايع‎ > saturation vapor pressure. 


This describes the behavior of e¿ asa Tumetien of T, 
assuming water vapor to be an ideal gas. It cannot be 
integrated exactly to give e, asa fūünctťtilion of T, since 
L(T) is not known to sufficient accuracy at more than a 
few temperatures [Weinreb, 1971]. 

The Goff/Gratch formula (Eq. 2) is an approximate 
solution of Eq. (1) considering the deviations from a 


perfect gas based on modern experimental data [List, 1963]. 


109,9 es = -7.90298 (T./T-1) 21209202908 log, y (T,¿/T) (2) 
Eee oo lie) 1) 
o (op 219 (.17-D . y 
+ log, ۾‎ Es 
where 
Lr. = steam point temperature (373.16°K) 
T = absolute (thermodynamic) temperature (°K) 
= saturation vapor pressure over a plane surface 
of pure ordinary liquid water (mb) 
ee Saturation pressure of pure ordinary liquid 


water at steam point pressure (mb). 
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Two saturation vapor pressures were calculated for 
each grid point using (a) the analysis-model field, 
giving ESAIR, and (b) the prognostic-model field, giving 
ESX. Then relative humidity parameters were calculated 


as follows: 


_ EAIR . 
RHR = ESAIR 100 
EX 
and RHX = ESX 100. 


3. Beta Visibility Parameter 
The computation of this parameter starts with the 
meeauction of an extinction coefficient, B, which is a 


function of windspeed and relative humidity. 
B = F(VVWW):-F (RHR or RHX) 


where VVWW = surface windspeed (m/sec) and 
RHR or RHX = relative humidity, 
and 


F(x) = A, + x (A, + x (A, + x (A, + x (A, + 200٠١ 


If the relative humidity input has a value greater than 


99.5 then it is set equal to 99.5. 
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The coefficients are as follows: 


For VVWW < 7 m/sec 


VVWW RHR or RHX 
A, 0.8065629 -0.4072407 x 10" 
A, 0.4852030 x 10! 0.3865717 
A, 0.5359734 x 10 ^ 50 140556 x l0 ^ 
A, 0.0 0.2496362 x 10°” 
A, 0.0 N OE 
A. 0.0 0.7388672 x 10 ? 
For VVWW > 7 m/sec 

VVWW RHR or RHX 
A, -0.8504248 x NOE los 
A, o O 0.583962 
A, -0.6052896 lassen. or 
A, TO o MO 
A, -0.1915719 x 10°° -0.328404 x 10? 
A, 0.3078907 x 10 ^ 0.1120986 x 107” 


Next, a new extinction coefficient is computed as, 


Bror = 8 + S where S is given as follows 
S Present Weather Code 
0.0 <50 
0.35 50-59 
0.2 60,61,80 
0. 6 626,81 
119 64,65,82 
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The scheme does not apply if weather codes other than 
those listed above are observed. The weather codes are 
defined in the Federal Meteorological Handbook No. 2 

[U.S. Departments of Commerce, Defense, and Transportation, 
1969]. 


Next, beta visibility is computed by 





BVISR = — , using RHR, and 
TOT 
3% 

BVISX = LI pe Sing RIX: 
TOT 


The modified beta visibility for use with prognostic times 
is computed without the weather code input by using the 


formula 


MBVIS = = ; 


and here RHX only is used for the relative humidity input. 
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APPENDIX C 


STATISTICS 


l. The coefficient of part determination, m may be 
interpreted as the proportion of the variance of the 
predictand that is explained by the regression equation. 


The computation of R* follows [Hill, 1979]. 


Y: = observed value of the dependent variable for 
case i. 

A = regression-specified value for case i 

Y = mean of the dependent variable 

(Y, -Y,) = residual for case 1, also called forecast 
error 

) (Y, Y) * = sum of squares about the regression line 

1 

) (v, -Y)* > Sum O! SGUarEEeS Of deviations about the mean 

1 

R = correlation coefficient between Y: and Y. 

pe = proportion of the variance of Y; that is 


"explained" by using Yar or 
=. 2 و‎ 2 
I SEE Y 


=, 2 
IO) 
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2. The F-to-Enter criterion used to enter variables in 
the stepwise regression procedure is given as follows 
Borel, 1979]. 

For each independent variable, Ku that is not in 
the equation at step (j+l), (j variables have already 


entered the equation); 


SES Enter = 
) (residuals at step j) - ) (residuals at step (j+l) 
1 1 


with x, in the e 


k 
) (residuals at step (3+1) with X, in the Suede) / 


1 (n-j-2) 


n = number of cases 


The F-to-Enter statistic is generally a measure of 
the importance of one variable relative to another. 
3. The goal in regression is to find the line, y, such 
that the sum of the squared residuals [) (v, -1)*] is 
minimized [Hill, 19791. For the line to be useful, it 
is required that the deviations between the observations 
and the line be smaller than the deviations between the 


line and the overall mean. Therefore, the quantity 
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Y (Y, -Y)^] -) (Y,-Y,) 1 should be large or one could 
say a good line has Y (Y-y) * small compared to (Y, Y) * 


y 
Yi 
, 
A 
(Yi- V) 
(Yi- Y) 
e ($-V) 
y gus ERO. EMEN esi 
^ 
y 
bo 
X 
The regression line is Y - مط‎ + biX, or generally, 
HU : b.X.; 


4. When an independent variable has a low tolerance it 
should not be included in a regression equation because 
its value can be expressed fairly well using a linear 


combination of variables already entered in the equation. 
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A variable with a low tolerance does not add significantly 
to the accuracy of a regression equation and may cause 
munerrcal and statistical accuracy problems [Hill, 1979]. 


The tolerance is computed by 


E 2 
TOLERANCE - D > E 
k 
where R is the multiple correlation coefficient of the 


entering variable, X with the set of independent 


les 
variables already in the equation, %£. If the computed 
value of tolerance is less than a preselected limit 

value, a prospective predictor cannot be selected for 


the regression equation as it is too highly correlated 


with the predictors already selected. 
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APPENDIX D 


VERIFICATION SCORE FORMULAE 


1. The two scores, percent correct and Heidke skill 
score, use a verification matrix as follows: (A 2x2 
matrix is used as an example, but the technique may be 


applied to any size matrix.) 


estimated 





1 = A+B 
k = C+D 
3 j = A+C 
> & = B+D 
بم‎ 
0 
Y 
e 
O 
_ A+D 
(a) Percent Correct = 28 x 100 
number of correct estimates 
total number of estimates 
(b) Heidke skill score = (A+D) - EXP 


(A+B+C+D) - EXP 


number of correct estimates - 


correct number expected due to chance 
total number of estimates - 


correct number expected due to chance 


_ (ij) + (2) 
م‎ 24 87+ 


ug 





E Bias Calculation 
Bias in estimating a given category = 


number of estimates of a given category 
number of observations of same category 


such as a or A 1 
af k 
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APPENDIX E 


SELECTED VERIFICATION MATRICES 


The following verification matrices show the number 
of observations in relation to the number of regression 
estimates for each visibility category, The top number 
ın each block is derived from dependent data and the 
bottom number from independent data. Row and column 


totals are given in the margins. 


l. Verification Matrix for 5P00: 


Regression estimated category 


i 01 ١ IV V 
2 2 174 2 70 521 
8 2 25 293 80 608 
9 4 5 133 zsnl 74 447 
i 2 1 99 165 60 327 
5 3 2 110 EE 150 588 
> 1 2 105 239 197 544 
2 1 0 58 299 698 
0 1 48 234 691 





e لن‎ 
O | حبر‎ 
o O 


0 54 455 16 1825 
1 39 448 2009 2557 


10 9 529 1581 1950 
TOTALS 
12 6 516 1379 2519 
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PEE Verification Matrix for 5P24: 





Regression estimated category 





I II Ia IV V 
21 21 IN 331 57 541 
I 
13 4 129 337 97 580 
> 13 10 91 269 81 464 
= 
O TI 
3 6 5 58 174 68 31 
لد‎ 
E 14 64 305 201 588 
EU III 
0 3 54 231 226 520 
ك‎ 
Q 2 34 260 398 698 
EN IV 
0 0 33 198 436 671 
3 0 B 410 1360 1804 
V 
3 3 44 350 2088 2488 
53 39 331 1575 2097 
TOTALS 
25 22 318 1290 2915 
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3. Verification Matrix for 5P48: 


Regression estimated category 


I X TIT IV V 
2 2 70 365 27 566 
0 0 40 336 n En 
1 1 2 34 286 143 466 
M 
8 0 0 24 147 239 302 
0 
a 1 33 | 269 295 598 
E 2 14 193 295 504 
z 
F 206 
u 
de 
O 


18 461 685 

IV 
9 175 468 652 
17 298 1472 L707 

V 
14 276 2126 2416 


3 5 172 1424 2498 





TOTALS 
0 2 101 1127 3191 





4. Verification Matrix; Probabilistic vs. Categorical 
This verification matrix shows results from dependent 
data for the probabilistic scheme of Aldinger [1979] vs. 
the SCAT categorical scheme of this study. The upper 
values in each block are for the probabilistic scheme, 


the lower values are for the categorical scheme. 


| Regression estimated category 





T TT 1 IV V 
9 


EI 

106 275 13 25105. 714 

i 7 3 106 504 94 714 
6 


00 
H 


5 76 275 264 198 93 906 
H II 
o 5 2 100 644 155 906 
لد‎ 
S 83 284 483 461 0649 1452 
TET 
2 2 90 902 456 1452 
> 
0 2 232 380 976 246 1911 
2 IV 
Q 
O 1 1 60 820 1029 E 
117 2 Bas 2240 1120 4137 
Y 0 1 53 1110 2973 4137 
459 | 1393 1599 2988 1681 
TOTALS 
15 9 409 3980 4707 
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